Decoherence through Spin Chains: Toy Model 
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The description of the dynamics of closed quantum systems, governed by the Schrodinger equation 
at first sight seems incompatible with the Lindblad equation describing open ones. By analyzing 
closed dynamics of a spin-i chain we reconstruct exponential decays characteristic for the latter 
model. We identify all necessary ingredients to efhciently model this behavior, such as an infinitely 
large environment and the coupling to the system weak in comparison to the internal couplings in 
the bath. 
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Quantum mechanics is distinguished from other prob- 
abilistic theories by the description of the relations be- 
tween various observables. These relations are described 
by coherences, off-diagonal entries of a matrix represent- 
ing a state. This fact has important consequences for 
dynamics of quantum systems. 

Understanding this dynamics has given rise to many 
important applications, such as parametric down- 
conversion m, perfect state transfer or sudden death 
Q (rebirth [J|, etc.) of entanglement. However this evo- 
lution might be seen at many different levels. The first 
is statical. One just applies a completely positive map to 
an initial state to get the final one [5-7]. This is a black 
box approach, where we are not interested in underly- 
ing physical phenomena. The second level is to solve the 
Schrodinger equation ^] to be able to trace the evolution 
of the system at any given instance of time. We assume 
here, however, that a quantum system is closed, which is 
often clearly false in light of experimental data. We can, 
of course, incorporate the environment, the "invisible" 
part of the system into the equation, which most often 
greatly comphcates its solutions. 

The third description of dynamics is based on the mas- 
ter equation Q. It has terms describing the free evo- 
lution governed by the Schrodinger equation and some 
phenomenological terms describing the interaction with 
the environment. 

While the completely positive map approach is gen- 
eral, the Schrodinger equation and the master equation 
seem incompatible. The former generates a unitary evo- 
lution group. This causes the probability of a survival 
of the initial state to diminish quadratically with time 
(for short times). Consequently, we observe the so-called 
quantum Zeno effect [lol . . If we frequently measure in 
the basis containing the initial state, the free dynamics 
of the system freezes - the survival probability of the ini- 
tial state tends to I in the limit of continuously repeated 
measurements. Moreover, while the free evolution is re- 
versible and even in time, the master equation describes 
exponential decays. The reverse evolution of an open 
system might make the state unphysical. How can these 
descriptions be agreed with each other? 

Here we want to show, how exponential decays of co- 
herence might follow from the Schrodinger equation. The 




FIG. 1. (Color Online) In an imaginary system, whose evolu- 
tion is governed by only two mutually irrational frequencies, 
1, TT, the survival probability Pi = + cos 2t + cos 2Trt) (left 
panel) never reaches 1 for t > and has a peak above 0.9 much 
earlier than P2 = ^ (3 -I- cos 2t + cos 2Tvt + cos 2et) correspond- 
ing to a system with three mutually irrational frequencies, 1, tt, 
and e (right panel). Only a system with infinitely many such 
frequencies allow (but does not guaranty) that such peaks 
never occur. 



evidence will be provided by a simple system, the semi- 
infinite chain of qubits coupled by the xx interaction. 
We will utilize the Heisenberg picture formalism [l2| - [l3 | , 
which turned out to be successful in analyzing the state 



initiation- free perfect state transfer (isl Ip 



and Hamil- 



18| 



tonian tomography with limited access 

B 

There have been already attempts to simulate deco- 
herence with a collection of spins, for example [l^. It 
should be pointed out, however, that Wang, Wang, and 
Su, as well as many other Authors 124- 34] , simulate the 
environment with a finite number of spins. In systems 
of finite size there is no possibility to have no remains 
of the initial state after infinitely long times. One may 
even expect a revival of the initial state with an arbitrary 
fidelity due to Poincare recurrences [20|, l2l| , however rare 
(see Figure 1). The deficits are removed in our approach. 

A detailed analysis of infinite spin chain baths, though 
in a slightly different context, in Refs. [13, Hj]. Ref. 
(2^ reconstructs the exponential decay in the approxi- 
mate calculations, while Ref. [2^ studies the behavior of 
the magnetization of one of the sites. These results are 
complementary to ours, which obtained within a different 
formalism. 

The system of our interest will be just a single qubit, 
labeled with 0. The role of the environment will be played 
by a wire of spins- labeled I to 00. Denote Pauli ma- 
trices acting on nth qubit by Xn , Yn , and Z„ . 
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are described by 



TABLE I. Number of n-step walks on non-negative numbers 
starting and ending at 0, and going through the 0th position 
k times. 



The overall Hamiltonian is 

JJ _ jjfree _j_ jjplug _^ jjwire 

= + ^{XoXi + YoYi) + y f2{Xa^+l + Y,Y,+i), 

(1) 

where we have taken the system of interest to be free, and 
given name "plug" to the interaction between the system 
and the environment (in analogy to "wire" , as the chain 
is called). We can now take the initial state of the 0th 
spin to be p = ^{1 + v- {Xg, Yq, Zq)), and the maximally 
mixed state as the state of the wire. In such a case, the 
evolution oi A = Xq , Yq , Zq and is given by 



A{t) 



E 

J=0 



(2) 



with L{-) — [H, ■]. In particular, we have 
Xo{t) + ao{t)Xo + ai{t)ZoYo + a2(t)^o^i^2 + ... (3) 

(confirm Ref. ^5|). Alternatively, we can write the evo- 
lution in the matrix form: 



(t) = cxp 
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(4) 



The problem can be hence mapped onto a classical ran- 
dom walk on the set of non-negative numbers. In our 
problem we observe only the 0th spin, hence we are inter- 
ested in the number of n-step (n being even) walks which 
start and end at 0. Since in general Jq ^ J, we need to 
categorize these walks in terms of how many times they 
pass the initial position (not including the start and the 
end, we call this number k). The numbers of such walks 
l{n, k) are given in Table 1. Collectively, these quantities 



l{n, k) 



(fc + l)(n-fc-2)! 
(f-fc-l)!f! 



(5) 



for 2k < n. Thus the original operators in the Heisenberg 
picture satisfy 

ao(i) = triXo{t)Xom = tr(ro(i)i"o(0)) 



2j 



2fi(l-j,2,2-2j,^) \ 
tr{Zo{t)Zom=tr{Xo{t)Xo{0)f. 



(6) 



Notice that the state of infinitively many qubits cannot 
be formally written as it faces the normalization problem. 
We have conveniently avoided this difficulty by studying 
the auto-fidelity of operators, ~1 < tr{A{t)A{0)) < 1. 

The remaining part of the paper will be devoted to 
analyzing how a general state of a qubit, p = ^{1+v^Xq+ 
VyYo + VzZq) behaves under the evolution of the whole 
chain. 

The first, trivial, case is of if = 0. Then the "environ- 
ment" contains only one qubit. We have ao(i) = cosKot. 
As expected, information contained in the system revives 
periodically, since it has nowhere to go. 

For K ^ 0, there are very few analytical solutions. 
One case is A' = where we have 
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Note hat the limit value of the Bessel function for its 
argument going to oo is 0. Hence, in the limit of long 
time the is no information at the 0th site about its initial 
state. It can be even argued that this information can be 
found nowhere specifically in the chain - it is completely 
delocalized. We already see how we benefit from consid- 
ering an infinitively high-dimensional environment. 

It is known from the theory or spin chains that the 
situation corresponds to an infinite spin 
chain, with the index ranging from —oo to oo. In such 
a case it is possible to analyze other phenomena. For 
example. Figure 2 presents the length of the Bloch vec- 
tor of the zeroth qubit, when initially it was in state 
|+)(+| = i(l + Xo), wheres the environment was fully 
magnetized (each spin of it in state |0)(0| = \{l + Zi)). 
Notice that the system periodically becomes completely 
polarized. 

Note that the relevant experimental results presented 
in (36j are in high agreement with our analysis, even 
though the number of parallel waveguides was only 21 
in the experiment. 



3 




FIG. 2. (Color Online) Evolution of the length of the Bloch 
vector v'^ = v'^ + Vy + of the 0th spin with Ko = V^K the 
initial state of the chain being | + 000...). v'^ = 1 corresponds 
to a completely magnetized state. 



The other analytical solution is for K 
we have 



Kq, for which 



ao{t) 



1 



{Ktf {Ktf {Ktf 
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144 



Ji{2Kt) 
Kt ' 



(8) 



Notice that this function also vanishes for t — )• oo, but its 
envelope is tighter than in the previous case. Instead of 
it now behaves like for long times. It is a general 
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the decay of ao will be 
= Kq will be referred to in 



trend: 

faster. The solution for K 
the later considerations. 

Further, not much can be said analytically. One thing 
we can show is that the first inflection point as we fix 
the characteristic time of the decay, , but go with , 
which characterizes the decay of information contained in 
the 0th spin, to infinity, we can show that the first inflec- 
tion point tends of ao{t) to 0. Namely, if one truncates 
d^aoi^x) to the second order 
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—^x^ = 0, 



(9) 



which has solutions at xo = ±-\/2 + ^) ^ j clearly 

tending to in the limit. This suggests that the quadratic 
approximation is relevant for shorter and shorter times. 
Figure 3 shows the following quantity: 



X 



{a'^ix) 



(10) 



in function of -S- 

^1 



ol'q{x) is Q!o(t) truncated up to 20th 
power of X. From this plot and Figure 4, which shows 
Q!o(t) for various ratios we conclude that in the 
limit of very high K the decay of information in the 
0th system becomes exponential. We must fail, however, 
to demonstrate it analytically. Since the Schrodinger 
equation gives only evolution even in time, the limit 

is exp ^— ^ , and the function has a point of non- 
smoothness at i = 0. We have ao{t) only as a Maclau- 
rin series, and hence we necessarily end up with a non- 
analytical function. 




FIG. 3. (Color Online) Logarithm of x given by Eq. (fTO)) in 



function of 




FIG. 4. (Color Online) ao(t) in function of time. The curves 
from the highest to the lowest aXt = 0.2 are plotted for = 

\/2, \/3, 2, \/5, \/6, \/7, 2\/2 and 2\/3. ^ has been fixed to 1. 



We can only describe the decay in a phenomenolog- 
ical way. Let us use the Trotter decomposition to the 
evolution. 



exp{iHt) = lim ( exp ( i!!£i j exp ' ^'"^ 



X exp 



(11) 



We start with the state of the 0th qubit p = |(1 -I- ^o)- 
The plug Hamiltonian changes operator Xq into ZqYi in 
some part. As we switch off the plug interaction and 
switch on the wire interaction, it instantly delocalizes 
ZqYi, practically mapping it to 0. In this approach we 
are unable to derive the lifetime of information, and we 
again face the problem of the Zeno effect. However, the 
obtained decays are clearly exponential. This allows us 
to speculate about more general models of decoherence. 

The specific Hamiltonian of the "wire" is in fact irrel- 
evant. The only thing that matters is that it maps any 
localized operator to an unlocalized one in an infinitely 
short time. Hence, two features of the environment are 
relevant. First that they are infinitely dimensional, and 
that the coupling constants between its constituents are 
scale properly with respect to the "plug" coupling con- 
stants. 

The choice of a specific plug determines the kind of 
decoherence. The main question is whether there exists 
an orbit of operators acting on the system, which com- 
mutes with the plug. For example, we could have taken 
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the Ising interaction, Hpiug = ZqZi^ rather than the xx 
interaction, we see that the rcsuhing map preserves the 
z component of the Bloch vector for the 0th spin (unless 
there is a transverse magnetic field acting on the sys- 
tem) . In this way we easily see that the "pancake map" , 
which maps a Bloch vector {vx,Vy,Vz) to {vx,Vy,0), is 
unphysical. For this observation, we do not need to re- 
fer to the concept of the complete positivity, but only 
to the non-existence of a plug, which commutes with Xq 
and Iq, but not with Zq. If from operators commuting 
with the plug we can build a Lie algebra, there exists a 
decoherence-free subspace. 

Of course, the toy model presented here allows to in- 
vestigate more phenomena related to decoherence. If our 
system consists of two qubits, each xx-coupled to a sep- 
arate wire, in the limit of — >■ oo, and we initialize the 
system in the singlet state, we observe that entanglement 
dies out after a finite time Q. If, however, we choose Kq 
sufficiently high in relation to K, ao{t) becomes oscillat- 
ing. We still may have the sum of squares of correlations 
between the two qubits, X^i j=i ^ij greater than 1, indi- 
cating entanglement [13] ■ We hence observe the so-called 
rebirth of entanglement occurring in non-Markovian 
environments, such as high-Q optical cavities. 

In conclusion, we have argued that a semi-infinite chain 
of spins- ^ can serve as a proper toy-model for decoher- 
ence. We have given the numerical evidence that in 
the properly taken limit of the weak system-environment 
coupling, exponential decays of information are observed. 
Moreover, due to a well established formalism of dynam- 
ics of spin-i chains, we have been also able to study some 
non-Markovian cases. We have captured all the neces- 
sary ingredients for decoherence (the environment of an 



infinite dimension and of an infinitely fast internal dy- 
namics) and discussed phenomena, such as the sudden 
death and rebirth of entanglement, or decoherence-free 
subspace. 

An important issue is whether the state of of the "wire" 
changes during the evolution. Of course, the unitary evo- 
lution preserves the amount of information contained in 
the whole system. As it leaks out from the 0th qubit, it 
appears in the "wire" . However, it would take a simulta- 
neous observation of infinitively many systems to recover 
this information, hence we can take the state of the wire 
effectively as maximally mixed. 

Interestingly, in the weak coupling limit we observe 
finite decay times for infinite values of Kq and K. 

Our result provokes some open questions. For exam- 
ple, we have conveniently assumed that spins in the en- 
vironment are already in a maximally mixed states. It 
shall not be difficult to perform a similar analysis for 
a completely magnetized environment, but can one effi- 
ciently model finite temperatures by biased random ori- 
enting spins? This would allow to generate mixed states 
of the system and the environment, while keeping the 
"universe" in a pure state. Then again, there would be 
a need to explain mechanisms behind this random orien- 
tation. Another question is whether there are any dif- 
ferences when a Heisenberg chain of spins- 1, not an xx 
chain of spins- ^, is used. The former kind has a gap be- 
tween the ground state and the first excited state, even 
in the thermodynamical limit. 
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